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ABSTRACT 

We report the discovery of a trend of increasing barium abundance with decreasing age for a large 
sample of Galactic open clusters. The observed pattern of [Ba/Fe] vs. age can be reproduced with 
a Galactic chemical evolution model only assuming a higher Ba yield from the s-process in low-mass 
stars than the average one suggested by parametrized models of neutron-capture nucleosynthesis. We 
show that this is possible in a scenario where the efhciency of the extra-mixing processes producing the 
neutron source ^'^C is anti-correlated with the initial mass, with a larger efficiency for lower masses. 
This is similar to the known trend of extended mixing episodes acting in H-rich layers and might 
suggest a common physical mechanism. 

Subject headings: Galaxy: disk, abundances, evolution; open clusters and associations: general, abun- 
dances 



1. INTRODUCTION 

Barium is one of the so-called heavy neutron-capture 
elements. It is formed by neutron captures occurring 
through two major mechanisms: (i) the slow capture 
process (s-process) and [ii) the rapid one (r-process). 
The r-process is thought to occur in high-mass stars, 
M > 8 Mq, exploding as type-II SNe (e.g., Kratz et al. 
2007; Travaglio et al. 1999). The s-process instead comes 
from low- and intermediate-mass stars, 1 M© ^ M < 
I 3 Mq (e.g., Busso et al. 2001), thanks to the repeated 
neutron expositions guaranteed by the activation of the 
neutron sources ^^C(Q;,n)^^0 and ^^Ne(Q;,n)^^Mg. The s- 
process contribution from low-mass stars had a dominant 
role in the solar mixture of Ba nuclei, so that in the 
Sun Ba comes from the s-process at about the 80% level 
(Arlandini et al. 1999). 

Barium abundances have been measured in stars with 
different metallicities (e.g., Burris et al. 2000 for halo 
stars; Cescutti et al. 2006 and Bensby et al. 2005 for 
thin and thick disk stars), constraining the nucleosyn- 
thetic processes at the origin of Ba. However, several 
details of the Ba production are still a matter of debate 
due to the free parameters affecting the models (see, e.g., 
for the r-process Mathews et al. 1992; Pagel & Taut- 
vaisiene 1997; Travaglio et al. 1999; for the s-process 
Busso et al. 2001). In the s-process, the main source 
of uncertainty is the still unknown physical mechanism 
promoting the penetration of protons into the He inter- 
shell zone at dredge-up; this penetration is necessary for 
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explaining the formation of the main neutron source ""^ C. 

In this Letter we address the evolution of Ba with age 
by means of Ba measurements in a sample of open clus- 
ters (OCs) and by comparing the empirical pattern with 
the results of a Galactic chemical evolution model. We 
find an unexpected trend, in apparent confiict with the 
currently adopted average Ba yields. We also identify 
in which direction the average yields should be corrected 
and why this is physically justified. 

2. SAMPLE AND ANALYSIS 

Our sample (see Table 1) consists of 20 OCs spanning a 
wide interval in age (~30 Myr-8 Gyr), Galactocentric ra- 
dius (7-22 kpc), and metallicity (-0.3 < [Fe/H] < +0.4). 
We analyzed a large collection of high-resolution spectra 
obtained by our group during the years with different 
instruments: UVES in sht-mode and FLAMES/UVES 
on the ESO VLT, CASPEC on the ESO 3.6m telescope, 
the Red-Long camera echelle spectrograph at CTIO, and 
SARG on the Italian National Telescope Galileo. The 
spectra of unevolved members were analyzed in the first 
10 OCs in Table 1, while the analysis of the older and 
more distant clusters is based on giant stars. 

Barium abundances were computed by means of equiv- 
alent width (EW) analysis and by using the driver blends 
in MOOG (Sneden, 1973 -2002 version). The grid of 1-D 
LTE model atmospheres of Kurucz (1993) was employed. 
Our analysis uses two Ba lines: AA 5853 and 6496 A. As 
well known, Ba lines are affected by hyperfine-splitting 
(hfs); thus we included in our line list the hyperfine 
structure and isotopic splitting, following the approach 
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by McWilliam (1998) and using log gf values published 
there. An 81% s-process composition was adopted for the 
Sun and the sample stars. Stark and radiative broaden- 
ing were treated in the standard way and for coUisional 
damping we adopted the classical Unsold (1955) approx- 
imation. 

As a first step, we derived the solar Ba abundance 
using the same method as for the sample stars: for 
the Sun we adopted Toff = 5770 K, log g=4.44 and 
^=1.1 kms/s (Randich et al. 2006), obtaining log n(Ba)= 
2.22 and logn(Ba)= 2.26 for A = 5853 and 6496 A, re- 
spectively. The average solar value log n(Ba)= 2.24 ± 
0.02 agrees very well with the meteoritic abundance of 
logn(Ba)=2.22 (Grevesse et al. 1996). Our analysis is 
strictly differential with respect to the Sun: we ap- 
plied a line by line solar correction, canceling out the 
uncertainties on the adopted set of atomic parameters 
(log 5/). The feature at A = 6496 A always yields higher 
ab undan ces, suggesting, as also stressed by Prochaska et 
al. (j2fl0l . the presence of an unidentified blend; however, 
the differential nature of our analysis should minimize 
this effect, in particular for dwarf stars with close-to- 
solar stellar parameters. Ba abundances were obtained 
adopting stellar parameters and uncertainties derived in 
our previous studies. 

Errors in Ba abundances for each star are due to un- 
certainties in EWs and stellar parameters. The stan- 
dard deviation from the mean of each cluster instead 
provides a good estimate of the global errors on the av- 
erage abundances"'^ . In order to estimate systematic er- 
rors, we compared our results with measurements from 
the literature for six clusters. We found a good agree- 
ment within the uncertainties (A[Ba/Fe]<0.06dex) for 
M67, Berkeley 20, and NGC 6253; the differences are 
instead larger for the Hyades (A[Ba/Fe]= —0.2 dex), 
Collinder 261 (A[Ba/Fe]=-|-0.2 dex), and Berkeley 29 
(A[Ba/Fe]= +0.1A dex). These results show that, 
whereas discrepancies with the literature are present, 
these are not systematic; in particular our Ba abun- 
dances are not systematically higher. Also, given the 
old age of Cr 261 and Be 29, our higher [Ba/Fe] ratio for 
these clusters does not affect our final conclusions. 

3. RESULTS 

The results are summarized in Table[Tl where we report 
the average [Ba/Fe] ratios for each cluster along with the 
standard deviation from the mean. In Figure 1 we show 
[Ba/Fe] ratios as a function of cluster age; different sym- 
bols indicate dwarf and giant members. In the figure we 
also show the solar [Ba/Fe] as well as our own measure- 
ment for a Centauri perfor med by using EWs publishe d 
by Porto de Mello et al. (jPorto de Mello et al. 20081 ). 
Our Ba abundance is in excellent agreement with this 
previous result (A [Ba/Fe] < 0.02 dex). The trend of in- 
creasing [Ba/Fe] from the oldest clusters to the youngest 
ones is evident. The figure also shows that [Ba/Fe] ratios 
are systematically higher for clusters where giants have 
been analyzed. We suggest that this might be due to the 
fact that our analysis is not strictly differential for gi- 
ant stars, as their parameters are different from the solar 

^ The line list with the EWs measurements as well as stellar pa- 
rameters, and Ba abundances of each single star along with errors 
are available upon request. 
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Fig. 1. — Average [Ba/Fe] as function of age for the sample 
clusters; filled and open squares denote clusters whose analysis is 
based on dwarf and giant stars, respectively. We also show the 
solar value (solar symbol) and that of the a Centauri (star with 
errorbars). 

ones. 

First, the question arises whether the observed trend of 
increasing [Ba/Fe] with decreasing age is real or rather 
due to any effect that, for a given intrinsic Ba abundance, 
would increase the observed EWs. In particular, mem- 
bers of the youngest clusters in our sample (ages below 
~ 500 Myr) are characterized by strong magnetic fields 
and are chromospherically very active; thus Ba II line 
strengths might be affected by both a "magnetic inten- 
sification mechanism" (Landi deglTnnocenti & Landolfi 
2004) and by non-LTE (NLTE) effects. In order to check 
whether magnetic intensification could be a problem, we 
measured the EWs of the weak Eu II 6437 A feature in 
the sample stars. Magnetic intensification would cause 
a similar (if not larger) effect for Eu, whose hyper-fine 
structure is notoriously more complex (the number of 
Zeeman sub-levels is higher). We found that, at variance 
with Ba, Eu EWs do not increase with the age, thus 
excluding that this process affects Ba EWs. 

On the other hand, NLTE effects might be important 
for Ba II features, in particular, for the A6496 A line and 
for stars with close-to-solar metallicity (Mashonkina et 
al. 1999). It is very likely that NLTE effects are ampli- 
fied in young stars due to the presence of hot chromo- 
spheres. While we cannot quantify this effect, we believe 
that [Ba/Fe] for clusters younger than ^ 500 Myr are 
likely overestimated. 

4. TIME EVOLUTION OF BA 

In Figure 2 we plot the average [Ba/Fe] of OCs (only 
those based o n dwa rf star analysis) and of field stars by 
Bensby et al. ( j2005l ) vs. stellar age. Although field stars 
show a larger scatter, the two patterns agree very well 
and our OCs do extend the same trend to younger ages. 
Qualitatively, the observed trend of Ba vs. age confirms 
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Fig. 2. — Average [Ba/Fe] as function of stellar age for the sub- 
sample of clusters whose analysis is based on dwarfs (filled circles 
and inverted triangles) compared with the abundance pattern of 
disk stars (open triangles) by Bensby et al. 1,2005. ). Filled trian- 
gles represent abundance measurements that probably need NLTE 
corrections. The model results are shown for two set of yields: 
(a) standard yields (Travaglio et al. 1999; Busso et al. 2001), long- 
dashed curve; and (6) enhanced s-process yields, dot-dashed curve. 
Both models show a peak at old ages due to the r-process from 
massive stars. 

the dominant role played, in Galactic disk stars, by the 
s-process occurring in low-mass red giants during their 
asymptotic giant branch (AGB) phase (Busso et al. 2001; 
Cescutti et al. 2006). For a more quantitative compari- 
son, we also show in Fig. 2 the time evolution of [Ba/Fe] 
computed with a Galactic chemical evolution model for 
the Milky Way. The model adopted is a gener alization of 
the multi- phase model by Ferrini et al. ()1992l ). originally 
built for the solar neighborhood, and later extended to 
model the entire Galaxy (Ferrini et al. 1994). We refer 
to Magrini et al. (2007, 2008) for a detailed descrip- 
tion of the model. In our "standard" model we adopted 
the average s-process and r-process yields by Busso et 
al. ()2001[ ) (extr apolat ed as in Cescutti et al. 2006) and 
Travaglio et al. (|2001h . respectively. While the r-process 
yields do not affect the [Ba/Fe] ratio in the metallicity 
range of our OCs (the r-process is important only for 
[Fe/H] < —2.2), the values of the s-process yields are 
crucial to reproduce the increasing trend of [Ba/Fe] with 
decreasing age of the cluster. 

The figure clearly shows that the model with stan- 
dard yields does not reproduce the empirical distribu- 
tion, and, in particular the observational evidence that 
at recent epochs Ba is sensitively enhanced both for field 
stars and for OCs. The standard s-process yields do not 
even fit the solar [Ba/Fe] value and result in a Ba abun- 
dance ratio [Ba/Fe] « —0.15 at 1 Gyr, to be compared 
with the measured value of ~ 0.2. In the original work 
from where the s-process yields were taken (Busso et al. 
2001) it was underlined that at every epoch and for ev- 
ery group of stars there was evidence of a large spread 
in their efficiencies. However, a simply randomly dis- 
tributed spread seems not to be the case: it would gen- 
erate from each given stellar generation a yield close to 
the average one adopted by us, and this is clearly inade- 
quate to explain the data. Instead, we "empirically" veri- 
fied that the observations, showing a very late increase of 



Ba in the Galaxy, would require an efficiency of produc- 
tion strongly weighted toward lower masses. Namely, the 
observations can be well fitted by increasing the yields 
of stellar masses between 1 and 1.5 Mq by a factor of 
~ 6 with respect to the extrapolation obtained from the 
average yields (roughly corresponding to the case ST/ 1.5 
in Busso et al. 2001). 

The results of a chemical evolution model where we 
have adopted "enhanced s-process" yields are shown in 
Fig. 2. In this model the solar Ba abundance can be re- 
produced, with contributions from the s- and r-process 
of 84% and 16%, respectively. As mentioned, the model 
with standard yields has instead a lower percentage of 
Ba produced by the s-process at the time of formation 
of the Sun. More in general, the model with enhanced 
yields provides a recent production of Ba in better agree- 
ment with the observations of OCs older than ^ 500 Myr. 
Both models fail to fit younger clusters and, indeed, it 
would be difficult to imagine a process capable of creat- 
ing Ba in the last 500 Myr of Galactic evolution, unless 
local enrichment is assumed. We suggest that the dif- 
ference between our model predictions and the measured 
[Ba/Fe] ratios for young clusters provides an estimate 
of the NLTE effects that affect the Ba determination in 
these stars. 

We finally stress that the classical plots of [Ba/Fe] 
against [Fe/H] (see, e.g., Travagho et al. 1999, 2001; 
Cescutti et al. 2006) suffer for a strong degeneration in 
terms of age, since they compress the most recent pe- 
riod. Conversely, the plot of [Ba/Fe] against age favors 
the youngest ages and thus is ideal to approach our prob- 
lem. We also investigated the presence of a possible re- 
lationship of [Ba/Fe] with metallicity among our sample, 
finding that no significant trend is present between these 
two quantities. 

5. DISCUSSION 

Our results show for the first time there is a need for an 
enhanced Ba production in low-mass stars, mainly from 
1 to 1.5 Mq, formed from ~ 10 to 5 Gyr ago (e.g. 
Charbonnel et al. 1996). We have shown that, when 
adopting the average yields that have been commonly 
used so far. Galactic evolutionary models are in sharp 
disagreement with data. An enhanced Ba production, 
like the one assumed in the previous section, seems to 
be necessary, but one has to verify whether this ad-hoc 
data-fitting is physically acceptable. 

The most promising possibility to obtain enhance- 
ments in the Ba production, especially in low masses, is 
that the intrinsic spread found normally in the s-process 
efficiency at any metallicity be not really a spread, but 
rather a dependence of the efficiency itself on the initial 
mass, in the sense of higher production for lower masses. 
As time passes, lower mass stars with larger Ba produc- 
tion start to contribute, and the average yield changes 
in time. Specifically, the formation of the neutron source 
^■^C requires a penetration of protons below the formal 
convective border at dredge-up. This is a typical extra- 
mixing process, like those necessary to account for the 
carbon and oxygen isotopic mix in low-mass red giants, 
and for the production and destruction of Li in them. 
These mechanisms are known to be increasingly more 
efficient for decreasing stellar mass. The evidence com- 
ing from the Ba enhancement indirectly suggests that 
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the mixing processes driving the formation of the neu- 
tron source behave similarly, and are maybe originated 
by the same cause (e.g. Busso et al. 2007). In general, 
we verified that the efficiency in s-processing required in 
very low mass stars in order to explain the data is close 
to the upper limit (case ST*2) suggested by Busso et al 
(2001), not to the average of the spread there considered 
(case ST/1.5). This must come together with the (al- 
ready known) increase in the dredge-up efficiency at low 
metallicity, so that this large amount of extra Ba can be 
carried to the surface and ejected into the ISM. 

Other possibilities were also considered such as: i) 
an increase in the initial s-element abundances of the 
contributing stars, due to previous enrichment by mass 
transfer from more massive companions; and ii) an in- 
crease in the r-process seeds on which slow n-captures oc- 
cur, assuming that old disk stars were still characterized 
by high r-process abundances. Both attempts proved 
incapable of increasing substantially the Ba yields. 

6. CONCLUSIONS 

High-resolution observations of OCs allowed us to de- 
rive Ba abundances in a uniform way in systems covering 
a wide age range. The [Ba/Fe] derived in OCs (but also 
in field stars, by other authors) shows a steady increase 
with time. In the framework of Galactic chemical evo- 
lution, this can be explained by an enhanced s-process 
production of Ba in very old AGB stars below 1.5 Mq. 

This need can be reconciled with the present under- 
standing of the s-process if an anti-correlation exists be- 
tween the initial mass of the star and the effectiveness 
of the extra-mixing mechanism responsible for forming 
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the neutron source. Stars at the lowest mass limit for 
contributing would in this case have a higher efficiency 
in s-processing. Accounting for the observations then 
requires that this efficiency be close to the upper limit 
indicated by Busso et al. (2001), i.e. an amount of ^'^C 
burnt per interpulse higher by a factor of three than the 
gross average (case ST/1.5) discussed in that work: this 
would produce an increase in the Ba yield by a factor of 
5-6. This scenario is plausible if the mechanisms driving 
the formation of the ^'^C pocket are similar to those pro- 
moting the slow circulations known to occur below the 
convective envelope. 

The youngest clusters of our sample, with ages < 
500 Myr, have even higher [Ba/Fe] ratios. We suggest 
that this is due to NLTE effects, and, in particular, to 
the effects of a strong, warm flux from a chromosphere. 
The comparison with our model suggests that NLTE ef- 
fects in young stars might increase the [Ba/Fe] ratios by a 
factor of two of more. Independent estimates of these ef- 
fects from NLTE calculations as well as a more complete 
analysis of the hypothesis here advanced of an anticor- 
relation between initial mass and s-process efficiency are 
strongly needed. 
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TABLE 1 

OC SAMPLE 



Cluster 


Age (Gyr) 


[Fe/H] 


Ref. 


[Ba/teJ 


# stars 


TeH (K) 


IC 2602 


0.035 


± 0.01 


1 


0.64 ± 0.07 


8 


4770-5760 


IC 2391 


0.053 


-0.01 ± 0.02 


1 


0.68 ± 0.07 


6 


4680-5970 


NGC 2516 


0.110 


0.01 ± 0.07^ 


2 


0.41 ± 0.04 


4 


5110-5659 


NGC 6475 


0.22 


0.14 ± 0.06 


3 


0.33 ± 0.11 


19 


5048-5888 


Hyades 


0.7 


0.13 ±0.05'' 


2 


0.30 ± 0.05 


4 


5079-5339 


Pracscpc 


0.7 


0.27 ±0.04 


4 


0.22 ±0.06 


6 


5720-6280 


NGC 3680 


1.4 


-0.04 ± 0.03 


4 


0.25 ±0.06 


2 


6010-6210 


IC 4651 


1.7 


0.12 ±0.05 


4 


0.20 ±0.03 


5 


5910-6320 


NGC 752 


2.0 


0.01 ±0.04 


5 


0.19 ±0.08 


9 


5531-6151 


M67 


4.3 


0.02 ±0.14 


6 


0.04 ± 0.05 


10 


5541-6223 



NGC 2324 


0.7 


-0.17 ±0.05 


7 


0.66 ±0.09 


7 


4300- 


-5100 


NGC 3960 


0.7 


0.02 ±0.04 


8 


0.50 ±0.06 


6 


4870- 


-5050 


NGC 2660 


0.8 


0.04 ±0.04 


8 


0.47 ±0.06 


5 


5030 


5200 


NGC 2477 


1.0 


0.07 ±0.00 


7 


0.40 ± 0.03 


6 


4950 


-5030 


Berkeley 20 


4.3 


-0.30 ±0.02 


10 


0.09 ±0.20 


2 


4400- 


-4770 


Berkeley 29 


4.3 


-0.31 ±0.03 


10 


0.40 ±0.10 


5 


4970- 


-5020 


NGC 6253 


4.5 


0.39 ± 0.08 


9 


0.14 ±0.10 


5 


4450- 


-6200 


Melotte 66 


5.4 


-0.33 ± 0.03 


10 


0.32 ±0.10 


6 


4717-4850 


Berkeley 32 


6.1 


-0.29 ± 0.04 


8 


0.24 ±0.15 


9 


4760- 


-4920 


CoUinder 261 


8.4 


0.13 ±0.05 


10 


0.22 ± 0.03 


7 


4350- 


-4720 



Note. — Results of [Ba/Fe] ratios for our sample OCs. We list; ages (Col. 2) from Magrini et al. (2008), with the exception of NGC 2516 
(WEBDA) and NGC 752 (Sestito et al. 2004); [Fe/H] values along with references for them and stellar parameters (Cols. 3 and 4); average 
[Ba/Fe] ratios along with the standard deviation (Col. 5); number of stars and range of effective temperatures (Cols. 6 and 7). Reference 
codes are as follows: l=D'Orazi & Randich (2008), 2=Randich et al. (2007), 3=Sestito et al. (2003), 4=Pace et al. (2008), 5=Sestito et 
al. (2004),6=Randich et al. (2006), 7=Bragaglia et al. (2008), 8=Sestito et al. (2006), 9=Sestito et al. (2007), 10=Sestito et al. (2008). 



